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We calculate the probability of electron-positron pair creation in vacuum in 3+1 
dimensions by an external electromagnetic field composed of a constant uniform 
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Pair creation of charged particles in vacuum by an external electric field was first studied by Schwinger 
'several decades ago [1]. Related problems have been discussed by many authors, for example, in Refs. 
[2-6]. Though half a century has passed since the publication of Schwinger's classic work, the subject of 
'pair creation remains a widely discussed one in string theory and black hole theory in these days. 
\ It seems that most authors of the cited works dealt with the problem only in electric fields. A magnetic 
f^*) .field, on the other hand, was left out in the cold. This may be due to the fact that a pure magnetic field 
i-H does not lead to creation of particle-antiparticle pairs. In spite of this fact, the presence of a magnetic 
OO .field does change the probability of pair creation by a pure electric field. This can be easily seen when 
'the magnetic field B is perpendicular to the electric field E. The situation may be changed by a Lorentz 
^ C| boost to another one with a pure electric field, provided |E| > |B|. The result is easily obtained, and 
ds obviously different from the case where B is absent. When B is not perpendicular to E, the problem 
O ,'is more complicated. The result for the general case cannot be easily obtained from previous ones and 
1^ .needs further study. As an exact and nonperturbative result can be achieved as in the case of a pure 
*^ 'electric field, the study is worthwhile and may be of some interest. 
J> . In this paper we consider a magnetic field B as well as an electric field E, both being constant and 
uniform, but with arbitrary magnitudes and directions. The field of electrons (or other charged fermions 



X 



. of spin i) is second quantized, while the electromagnetic field is treated classically as a background field 
^ for the electrons. The probability of electron-positron pair creation in vacuum is calculated exactly. More 
specifically, we first deal with the relatively simple case where B || E (B points at the same or opposite 
direction of E). For the general case, one can always find by making a Lorentz boost an inertial frame 
K' where the transformed fields satisfy B' || E'. In the system K' the result is obtained in terms of E' 
and B'. As the probability should be invariant under Lorentz boost, we obtain the result in the original 
system by using the relation between (E',B') and (E, B). When E = the probability vanishes while 
for B = it reduces to Schwinger's result, as expected. This is done in Sec. II. 

In previous works main attention was paid to the problem in ordinary 3+1 dimensions. In lower spatial 
dimensions the result might be expected to be somewhat different. In Sec. Ill we turn to the problem in 
2+1 dimensions. In this case the magnetic field has only one component and the electric field has two, 
thus the problem is simpler than in 3+1 dimensions. In Sec. IV we calculate the result in 1+1 dimensions. 
In this case there is no magnetic field and the electric field has only one component, so the problem is still 
simpler. The results in lower dimensions are similar to that in 3+1 dimensions, but cannot be trivially 
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obtained from the latter. The difference among them can be easily seen from the appearances of the 
corresponding results. These are summarized in Sec. V. 



II. 3+1 DIMENSIONS 



We use natural units where Ti — c = 1 throughout this paper. Consider the electron with mass m 
and charge e < (the following results are applicable to other charged fermions of spin i), moving in a 
background electromagnetic field described by the vector potential A^. The field of the electron is second 
quantized, while is treated classically. The vacuum- vacuum transition amplitude can be shown to be 
[1, 2] 



So = exp 



-Tr In 



7 • P — m + ie 
7 • (P — eA) — m + ie 



(1) 



where e = + , and the Tr indicates a complete diagonal summation over the space-time coordinates as 
well as the spinorial indices. In this section we consider the problem in ordinary 3+1 dimensions. But 
Eq. (1) holds in 2+1 and 1+1 dimensions as well. In Eq. (1) (the independent variables of A^) and 
P M are now operators satisfying 



Xn\x) = x^x), (xlP^ip) = id^(x\ip), 
where \ip) is an arbitrary state. Consequently 



(2) 



(3) 



where g^ v = diag(l, — 1, — 1, —1) and similarly in lower dimensions. Using the charge conjugation matrix 
and the fact that the trace of an operator is invariant under matrix transposition we have 



So = exp 



-Trln- 



7 • P + m — ie 



7 • (P — eA) + m — ie 
This holds in lower dimensions as well. Multiplying Eqs. (1) and (4) we have 



5$ = exp | -Trln 



ie 



[7 • (P - eA)} 2 -m 2 +ie. 
Taking the module we obtain the vacuum-vacuum transition probability in the form 



where 



w(x) 



Re < x 



\Sq\ 2 = exp 



trln 



— J dx w(x) 



[7 • (P - eA)} 2 - m 2 + ie 



ReW(x) 



(4) 



(5) 



(6) 



(7) 



is to be interpreted as the probability, per unit time and per unit volume, at the space-time position x, of 
electron-positron pair creation by the external electromagnetic field. For constant uniform electromagnetic 
field it is expected to be independent of x. In the above equation the tr indicates ordinary diagonal 
summation over spinorial indices. Using the identity 



, a + ie f ds . 
m = / — 1 

b + ie J Q s 



e is(fc+ie) _ e is(a+ie)-l 



and the ralation 



[7 • (P - eA)} 2 = (P- eA) 2 -a^F^ 



(8) 



(9) 



where 



2 



we have 



d^-d^, <7^ = -[ 7 ^,7l, (10) 



P Tfl -\- l€ t . ... _ /s[ n! - _ ( \)- I ut> j-h 1 /.s/'~ (11) 



-t - "I T _ f°° df_ _ is(m 2_. e) 

[ 7 • (P - eA)] 2 - m 2 + ie J s 



Up to this point wc have only reviewed some results obtained by previous authors [1,2]. These are 
necessary for further discussions. Wc emphasize that Eqs. (7) and (11) are valid in lower dimensions 
as well as in 3+1 dimensions. The following calculations depend on the spatial dimension, thus we 
will henceforth deal only with the (3+l)-dimensional case in this section. We will return to the lower- 
dimensional case in the following sections. 

Now we consider a constant uniform electromagnetic field where B || E. Without loss of generality we 
take 

E = Ee x , B = Be x , (12) 

where e x is the unit vector in the x 1 direction, E and B are constants which may be positive or negative. 
We have then for the nonvanishing components F i = E, F23 — —B, and 

-lia^F^ = 7 V £ - 7 2 7 3 B. 

Using the properties of the 7 matrices, we have 

(-Ita^F^) 2 = E 2 -B 2 + i2EB l5 , (13) 

where 75 = 7 s = i7°7 1 7 2 7 3 . As 7 | = 1, the eigenvalues of 75 are ±1, both being double degenerate. Thus 
the eigenvalues of (— ia^ F^/2) 2 are (E ± iB) 2 , both being double degenerate, and the four eigenvalues 
of -ia^F^/2 arc E±iB,-(E± iB). Therefore we have 

tr exp (-^iesa^F^) = 4 cosh(e£s) cos(ePs), (14) 

and 

The next step is to calculate the matrix elements (x\e ts< - p ~ eA ) 2 \x) and (x\e tsp2 \x). The second is easy. 
We have 

(x\e lsp2 \x) = J dk (x\e lsp2 \k)(k\x), (16) 

where 

= jt? e ~ ik ' x > m = ji^ elk ' x - (17) 

Using Eq. (2) it is easy to show that 

P li \k) = k li \k), (k\P^ = k^(k\. (18) 
Substituting Eqs. (17) and (18) into Eq. (16) we obtain 

The calculation of another matrix element is somewhat complicated. We denote 

X» = (X ,X,Y,Z), P" = (P ,P x ,P y ,P z ), (20) 
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and choose 

A" = (A Q , A x , A V ,A Z ) = (-EX, 0, 0, BY), (21) 
which corresponds to the field strengths (12), then 

(P - eA) 2 = (P + eEX) 2 - (P z - eBY) 2 - P 2 - P 2 . (22) 
It is not difficult to show that 

(P + eEX) 2 = e iPoP*/eE e 2 E 2 x 2 e -iP P x /eE^ (23a) 

(P z - eBY) 2 = e -iPyPJeB e 2 B 2 Y 2 e iP y P z /eB_ (23b) 

Substituting these ralations into Eq. (22) we have 

(P - eA) 2 = - e iP ° p */ eE (P 2 - e 2 E 2 X 2 )e-^ p *' eE - e-^ p ^ eB (P 2 + e 2 B 2 Y 2 )e^' eB , (24) 
and thus 

e is(P-eA) 2 _ e iP a P x /eE-iP y Pz/eB e -is(p3-e 2 E 2 X 2 ) e -is(p3+e 2 B 2 Y 2 ) e -iP„P x /eE+iP y Pz/eB ^5) 

With these preparations we write down 

(x\e ls ( p ~ eA ^ \x) = J dkdk' {x\k'){k'\e ts{p - eA ^\k){k\x). (26) 
Using Eqs. (25), (17), (18), and S[(k' x — k x )/eE] = \eE\S(k' x — k x ) etc., after some algebras we arrive at 

{x \ e is(P-eAf lx) = \^3^ tre ^s { P 2 -e 2 E 2 X 2 hYe -MP 2 + e 2 B 2 Y 2 ) ^ (2?) 



where 



tre -" (Pl-^E 2 x 2 ) = J dkx (k x \e-^ p *- e2E2x2 *>\k x ) 



and similarly for another trace, where k x is the first spatial component of fc M . By comparison with the 
harmonic oscillator or by using the technique of path integral, one can find the following results. 

tYe -MP 2 -e 2 E 2 x 2 ) = 1-—, (28a) 

2sinh(|e£|s)' V ; 



tre — (Pi+^B 2 Y 2 ) = J_^_ (28b) 

2sin(| eJ B|s) 1 ; 



Substituting these into Eq. (27) we obtain 



t x \ e *<P^A) 2 1 , i\eE\\eB\ 

[l l) ~ 167r 2 sinh(|ei?| S )sin(|ei?| S )- { ' 



Combining the results (7), (15), (19), and (29) we arrive at 



oo -tm s 



\eE\\eB\coth{\eE\s)cot(\eB\s) - \ 



(30) 



where the e-dependent term in the exponential has been dropped as it is no longer necessary. This is 
obviously independent of x as expected. Using the identity RelU = (W + W*)/2 and making the change 
of variable s — > — s in W we obtain 
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|e J B||eB|coth(|e J B|s)cot(|e J B|s) - 



(31) 



It is easy to see that the integrand has singularities (simple poles) at and ±nir/\eB\ (n = 1, 2, . . .) in the 
integration path. Thus the integral in the above equation is not well defined. An appropriate prescription 
should be employed. There are mainly three different prescriptions. The first is to replace the integration 
path by a straight line a bit above the real axis on the complex s plane, i.e., a straight line from — oo + ie 
to +oo + ie where e = + . The second is to replace the integration path by one from — oo — ie to +oo — ie. 
The third is to keep the original path but take the Cauchy principal value. It turns out that only the 
first prescription gives a physically acceptable result. The second or the third will give a negative result 
for w when B = 0, and thus are not acceptable. When B = the integrand has one singularity in the 
integration path, the simple pole at s = 0. An appropriate prescription is also needed in this case. It is 
the first prescription described above that was implicitly used in Rcf. [2] . The validity of the prescription 
was confirmed by the coincidence of the result with that of Schwingcr obtained by a different method. We 
adopt this prescription, and close the contour of integration at infinity by a semicircle in the upper half 
plane. The integrand has simple poles at nm/\eE\ (n = 1,2,.. .) in the upper half plane. The integral 
can be evaluated by using the residue theorem, and the result turns out to be 

e 2 \EB\ g 1 , / \B\\ ( nirm 2 

n=l 



-II - 2. - coth I mr— I exp 1-— I . (32) 



This is the result for the simple case B || E, as indicated by the subscript. The convergence of the series 
is obvious. It is easy the see that w — when E — 0, which means that a pure magnetic field cannot lead 
to pair creation as expected. For B = 0, on the other hand, we use cothw — > 1/u (u — ► 0), and obtain 



El / n7TTO 2 \ . . 

^ cx p "W ' (33) 



e E 2 \ - 1 / nirm 

n=l v ' 

where the subscript indicates a pure electric field. This is just Schwinger's result. When \B\ = \E\, we 
have 

e 2 E 2 1 , / nirm 2 \ 

w \b\=\e\ = 2^ - coth(nTr) exp I — r-jn ) ■ (34) 

n— 1 ^ ' ' ' 

Since coth(n7r) > 1 and 1/n > 1/n 2 , the series in Eq. (34) is obviously larger than that in Eq. (33), and 
the overall factor in Eq. (34) is also larger than the one in Eq. (33) by a factor n, hence w\b\=\e\ > kide- 
However, since \eE\ <C m 2 as currently available, the dominant term in both Eqs. (33) and (34) is the 
first term, and thus W|b|=|_e| ~ kwe- If B increases, the result may be still larger. However, If E cannot 
be significantly raised, the addition of a magnetic field cannot raise the probability greatly. 

We can now turn to the general case where both E and B have arbitrary magnitudes and directions. 
We can always find another frame of reference K' by a Lorentz boost where the transformed fields E' 
and B' satisfy B' || E'. In the system K' we can find the probability in terms of E' and B' by using Eq. 
(32). We know that the probability is a Lorentz scalar. This is because the total probability J dx w(x) is 
a Lorentz scalar as it only involves a process of number counting, and the space-time volume dx is also 
invariant under Lorentz boost. Hence the probability in the original system can be obtained by using 
the relation between (E, B) and (E',B'). Another approach to the result is also available. Since w(x) 
is a Lorentz scalar, it can only involve Lorentz invariants constructed from E and B. In the special case 
B || E, it must reduces to the result (32). This enables us to obtain the general result more easily. We 
define the Lorentz invariants 

T = E 2 B 2 , Q = 2E • B, (35a) 



then the result reads 
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e 2 £B^ 1 / B\ ( mrm 2 \ 
w = X, - coth [nn- j exp j . (36) 



4tt 2 

71=1 

Now we can examine the result by another special case where B _L E. In this case we have Q = 0, and 



£ = VE 2 B 2 , B = (37a) 

if |E| > IBI, or 



£ = 0, B = VB 2 - E 2 (37b) 

if |E| < |B|. The case (37a) is equivalent to that of a pure electric field and the result is given by Eq. 
(33) where \E\ is replaced by £ = \/E 2 -B 2 . The case (37b) is equivalent to that of a pure magnetic 
field and w = 0. These are all expected results and further confirm the result (36). 

III. 2+1 DIMENSIONS 

In this section we calculate the probability of pair creation in vacuum in 2+1 dimensions. This cannot 
be trivially obtained from the result in 3+1 dimensions. In two spatial dimensions the magnetic field 
has only one component, while the electric field has two. There is nothing like E • B. The only Lorentz 
invariant constructed from E and B is E 2 — B 2 . Consider an electromagnetic field (E, B) where both E 
and B are constant and uniform. There are two different situations to be distinguished. The first one is 
characterized by the inequality |E| > \B\, while the second by |E| < \B\. The first situation is equivalent 
to one with a pure electric field. One can calculate the result for the simple case with a pure electric field, 
and get the result for the general case by the method of Sec. II. This will be discussed in detail in the 
following. The second situation is equivalent to one with a pure magnetic field. Similar calculations give 
a vanishing probability in this case. This is an expected result and we will not discuss it in detail. 

In Sec. II we have emphasized that Eqs. (7) and (11) are valid in lower dimensions. We will begin 
with these equations. Consider a pure electric field E which is constant and uniform. Without loss of 
generality we choose 

E = Ee x , (38) 

where E is a constant which may be positive or negative. We have then for F^ v the nonvanishing 
component F i — E, and 

exp(-^ies(j^F^) = exp(e.Es7V )■ (39) 

As (7°7 1 ) 2 = 1, and tr(7°7 1 ) = 0, the trace of the above expression can be evaluated directly with the 
result 

tr exp (-\ies(j^ v F^) = 2cosh(e£s). (40) 

Note that in 2+1 dimensions the 7 matrices are 2x2 ones and thus trl — 2. Combining Eqs. (40) and 
(11) we have 

trln r + = 2 r -e-^ m2 -^[cosh(eEs)e^ p -^ 2 - e isP \ (41) 

[7 • (P - eA)} 2 -m 2 +ie J s 

The next step is to evaluate the matrix elements (x\e ls ^ p ~ eA ^ 2 \x) and (x\e lsp2 \x). The second one can 
be easily worked out with the result 

(x\e i ° p2 \x)= l~J 3 . (42) 
4(27r)2 S 2 

This is rather different from the corresponding result in 3+1 dimensions, but the calculation is similar. 
For the first one, we denote 
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X» = (X ,X,Y), P» = (P ,P x ,P y ), (43) 

and choose 

A" = (A ,A X , A y ) = (-EX, 0, 0), (44) 
which corresponds to the field strength (38), and results in 

(P-eAf = (P Q + eEXf -Pl-P*. (45) 

Using Eq. (23a) we have 

e is(P-eA) 2 = e iP P !c /eE e -is(Pi-e 2 E 2 X 2 ) e -iPoP !c /eE e -isP^ ^ 

The subsequent calculations are similar to those carried out in obtaining Eq. (27) but simpler, the result 
reads 

{x \ e HP-eA) 2 lx) = (l-«)je^| ^-eW)^ (47) 

V ' 17 2(2^)iVi V ' 



Using Eq. (28a) we obtain 



( X \e^ p -^ 2 \x) = j 1 "^ 1 . (48) 



4(27r)3 v ^sinh(|e£;|s) 
Combining Eqs. (7), (41), (42), and (48) we arrive at 



W(x) = — V / ds 



2(27r)2 Jo s'z 



|e£|coth(|e£|s) - i 



(49) 



If we are going to use the residue theorem for contour integrals, we must treat the integrand carefully 
since it is a multivalued function in the complex s plane. We cut the s plane from to —zoo along the 
imaginary axis, and define args = in the positive real axis, then the integrand is single valued in the 
cut plane. We use the identity ReW = (W + W*)/2 and making the change of variable s — > s' = e 47r s in 
W to yield 



1 + i [+°° 

w = o- / as 

4(27r)i J- 



4(27r)2 J-oo S2 



\eE\coth(\eE\s) - - 
s 



(50) 



Now the integrand has one singularity in the integration path, the origin s = 0. This is a branch point 
of the integrand. Moreover, the integrand tends to infinity like l/^/s when s — > 0. As the lower half 
plane has been cut, we have now two different prescriptions for the integral: to replace the integration 
path by one from — oo + ie to +oo + ie, or to take the Cauchy principal value. It turns out that the 
two prescriptions give the same result. We adopt the first prescription, close the contour of integration 
at infinity by a semicircle in the upper half plane and use the residue theorem to evaluate the integral. 
Note that there are simple poles mre™/ 2 /\eE\ (n — 1, 2, . . .) of the integrand in the upper half plane. The 
result turns out to be 

|e£ll ^ 1 / mrm 2 \ 

n—l v 1 1 ' 

where the subscript indicates a pure electric field. 

The case of a pure magnetic field B can be treated in a similar way. It turns out that the probability 
vanishes as expected. 

To conclude this section we consider the general case of an electromagnetic field (E, B) where E has 
an arbitrary direction. As pointed out at the beginning of this section, one must distinguish between the 
two different cases |E| > \B\ and |E| < \B\. The results can be easily obtained by the method of Sec. II. 
When |E| > \B\ we have 

|e|5£§ JS, 1 / mrm 2 \ . . 

n—l x ii/ 
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where 



£ = \/E 2 -B 2 . 



(52b) 



When |E| < \B\ we have a vanishing result, since the case is equivalent to one with a pure magnetic 
field. 



IV. 1+1 DIMENSIONS 

In 1+1 dimensions there is no magnetic field, and the electric field has only one component. Thus the 
problem is still simpler. As before, we begin with Eqs. (7) and (11). Consider the electric field 



E = Ee x , 



(53) 



where E is a constant which may be positive or negative. The nonvanishing component of F M „ is Fq\ = E, 
and 



cxp(-±ies<7^F^) = cxp(e£s7V). 



(54) 



This is exactly the same as Eq. (39). In 1+1 dimensions the 7 matrices are also 2x2 ones, thus Eq. (40) 
remains valid. With Eqs. (40) and (11), we have a result of the same form as Eq. (41) where P 2 , say, is 
of course different in different dimensions. It is easy to show that 



As before, we denote 



and choose 



(x\e isp2 \x) = 

x 1 1 ' 4tts 



X» = (X ,X), P» = (P ,P X ), 



A» = (A ,A x ) = (-EX,0), 
which corresponds to the field strength (53), and results in 

(P-eAf = (P + eEX) 2 -Pl 

Using Eq. (23a) we obtain 



It is now quite essy to show that 



\eE\ 



x 1 1 1 2n 



Using Eq. (28a) we have 



(x\e 



is(P-eA) 2 



\x) 



\eE\ 



47r sinh(|e^|s) ' 
Combining Eqs. (7), (41), (55), and (61) we arrive at 



1 r°° p — im 2 s 

W{x) = — I ds 



\eE\ coth{\eE\s) - i 



70 s 

Using the identity ReW — (W + W*)/2 and making the change of variable s 



— s in W we obtain 



1 



^ + 00 ^im s 

tr = — I ds 

4ir 



\eE\coih(\eE\s) 



1 



(55) 

(56) 
(57) 
(58) 
(59) 

(60) 
(61) 

(62) 

(63) 
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The integrand is regular everywhere in the integration path, thus no prescription is necessary here. We 
close the contour of integration at infinity by a semicircle in the upper half plane and evaluate the integral 
by the residue theorem. The result turns out to be 

\eE\ 1 / nirm 2 \ 

n— 1 v 1 17 

The convergence of the series is obvious. This can also be written in a closed form 

(64') 



V. SUMMARY AND DISCUSSIONS 

In this paper we calculate the probability of electron-positron pair creation in vacuum in external 
constant uniform electromagnetic fields in 3+1 and lower dimensions. The results are also applicable to 
other charged fermions of spin i. In 3+1 dimensions the most general result is given by Eq. (36), where 
E and B may have arbitrary magnitudes and directions. The result for the ralatively simple case where 
B || E is given by Eq. (32). For a pure magnetic field the probability vanishes, which means that pair 
creation cannot occur in this case, as expected physically. For a pure electric field, our result reduces to 
that of Schwinger, Eq. (33). In 2+1 dimensions the magnetic field has only one component while the 
eletric field has two. When |E| > \B\ the result is given by Eq. (52). For the special case of a pure electric 
field this reduces to Eq. (51). When |E| < \B\ the result vanishes since the situation is equivalent to one 
with a pure magnetic field. In 1+1 dimensions there is no magnetic field, and the electric field has only 
one component. The result is given by Eq. (64), or by Eq. (64') in closed form. 

For a pure electric field, the results in the several cases of different dimensions studied can be written 
in a unified form 

= (1 + **») {2n)d L ^TT)7I ex P ( " W J ' (65) 

n— 1 ^ 

where d = 1,2,3 is the spatial dimension, and Sd3 = 1 when d = 3 and vanishes in other cases. The 
additional factor 2 in 3+1 dimensions is due to the double spin states. More specifically, in 3+1 dimen- 
sions, with a given momentum there are two linearly independent solutions of positive energy and two of 
negative energy to the free Dirac equation, while in 2+1 or 1+1 dimensions there is only one solution for 
each energy. 

Pair creation of charged scalar particles in vacuum in external electromagnetic fields can also be studied 
in a similar way. In this case, however, the probability is reduced by the presence of a magnetic field, 
thus we do not discuss the problem in detail. 
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NOTE ADDED IN PROOF 

After this paper was accepted for publication, we became aware of some recent works dealing with the 
same or relevant problems by different methods [7-9]. We thank the authors for bringing their works to 
our attention. 
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